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A microscopic formalism is developed that includes the coupling to two particle-hole phonons 
in the particle-hole propagator by extending the dressed random phase approximation (DRPA) 
equation for a finite system. The resulting formalism is applied to study the low-lying excitation 
spectrum of ^^O. It is observed that the coupling to two-phonon states at low energy generates 
excited states with quantum numbers that cannot be obtained in the DRPA approach. Nevertheless, 
the two-phonon states mix weakly with particle-hole configurations and participate only partially 
in the formation of the lowest-lying positive-parity excited states. The stability of the present 
calculation is tested vs. the truncation of model space. It is demonstrated that when single-particle 
strength fragmentation is properly considered, the present formalism exhibits convergence with 
(N , respect to the chosen model space within the confines of the chosen approximation scheme. 
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; I. INTRODUCTION 

I Recent interest in the spectroscopic factors for the removal of low-lyingpyo andp^^ strength from ^^O points to 
a substantial discrepancy between experiment Q and theory [3 13. H. 13. ItII^I^ I1(1 fill . The contribution to 

the reduction of these spectroscopic factors due to short-range correlations is important and is uniformly calculated 
to be about 10% in ^^O using different realistic interactions and theoretical approaches. The best agreement with the 
O [ data is obtained when also long-range correlations are included in the calculations IE 1^ ■ The results of Ref. 
2 ■ demonstrate a particularly strong correlation between the theoretically calculated excitation spectrum of ^^O and the 
*~J " resulting one-hole spectral function, which describes the excitations of the residual A — 1 nucleus. The appearance 
^ . of an additional P3/2 fragment at small missing energy when a low- lying excited 0"*" state is obtained is one such 
k>( ] correlation 0. Another one is the appearance of low-lying si/2 and 1^5/2 strength when the first 3^ state in ^^O is 

. correctly described Both features are in accordance with the experimental data 1]. 

d The calculation of Ref. [^] employed a Random Phase Approximation (RPA) approach to describe the excitation 

spectrum of ^^O. The quality of the spectrum in this approach is, however, quite inadequate. Even with realistic 
G-matrices as residual particle-hole (ph) interaction one generates at most a few low-lying collective isoscalar states of 
negative parity (3~ and 1~) and no low-lying isoscalar positive parity states are obtained 0|. These considerations 
suggest that an adequate description of the sp fragmentation, at low energy, requires a better description of the excited 
states within the same framework of self-consistent Green's function (SCGF). 

This improvement of the SCGF method is also important because it can be applied to heavier nuclei. This is possible 
since the method describes the features of the low-energy spectrum in terms of interactions between a relatively small 
number of quasiparticles and collective excitations. These modes represent the experimentally observed excitations 
of the system and therefore they include both the effects of short- and long-range correlations. The fragmentation 
of the single-particle (sp) strength is thus self-consistently included. Early applications to ''^Ca and ^"^Zr have been 
reported in Refs. p^. Il6l|. These calculations did not include self-consistency but nevertheless were able to yield a 
reasonable description of giant resonances and Gamow- Teller states as well as the low-lying collective states. Moreover 
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a formalism based on the Bethe-Salpeter equation for ph excitations can be extended to include the effects of the 
continuum Although the SCGF method has only been applied in doubly closed-shell nuclei, it is feasible that 
increased computational power will allow applications to open-shell systems. An exploratory calculation for semi- 
magic isotopes of Ni and Sn was reported in 18]. In the near future a huge amount of data will become available for 
unstable nuclei from radioactive beam facilities. For this reason it is important to develop techniques that are flexible 
enough to describe all the above effects. 

In this paper, we report a first step toward an extension of the SCGF formalism in this direction. We chose to use 
^^O as a test system for the following reasons. First, a successful calculation for this nucleus will help resolving the 
aforementioned discrepancies for the spectroscopic factors. Second, this system has been studied by various approaches 
in the past. From the results of these calculations one can infer the relevant physical ingredients to be included in 
the SCGF method. Third, since the spectrum of "'^^O is rather complicated, meeting the challenge of describing this 
system will virtually guarantee good results when the same method is applied to heavier nuclei. 

Shell-model calculations for this nucleus have been reported in Refs. |19l I20I I21L which indicate the particular 
importance of four-particle four-hole (4p4h) admixtures to the OpOh ground state in generating the first excited 0+ 
state. Other positive parity states at low energy are dominated by 2p2h components, as indirectly confirmed by 
inelastic electron scattering | 2^ . Some of this resulting physics was anticipated in terms of deformation effects which 
simulate these types of many-particle many- hole admixtures 24]. A complementary point of view is given by the 
Interacting Boson Model of Refs. [25ll2^ . There, the low-lying positive parity states are understood in terms of the 
coupling of 3^ and 1^ isoscalar states. One should note that, due to the predominant Iplh nature of these excitations, 
the coupling of different phonons generates the npnh configurations relevant to this problem. In Refs. [25I |2^ the 
3~ and 1~ states were used as phenomenological boson degrees of freedom to generate the spectrum of the nucleus. 
The connection with the underlying fermionic description was indicated but not explored completely. This relation 
was considered in Ref. |27j |. There it was argued that the microscopic ph interaction contains two-phonon exchange 
contributions which include the actually observed low-lying states of "'^^O themselves, thereby generating the correct 
number of low-lying states observed at low excitation energy. 

In this work we begin to implement the physical ingredients proposed in Refs. '25', "2^ by extending the RPA to 
include the coupling to two-phonon states. The basic idea of the present approach is illustrated in Fig. ^ This figure 
depicts the coupling of a ph state to two intermediate phonons that are described by the ph propagator itself. If 
the two intermediate phonons have been computed using the RPA equations, they will already provide a reasonable 
description of the low-lying collective isoscalar 3~ and 1~ excited states. These phonons can be sufficient to generate 
the quantum numbers of the most important positive parity states. Particularly relevant is the coupling of two 3~ 
phonons, 3~ ® 3~ = 0+ © 2+ © 4+ © 6+, that represent some of the correlated 2p2h states of ^^O. In the framework of 
SCGF one employs dressed sp propagator in the construction of the microscopic phonons. This results in a Dressed 
RPA (DRPA) approach for the ph calculation and in its extension to the coupling to two-phonon states. In the 
present paper, we employ the self-consistent sp propagator of ^^O computed in Refs. 0,13. It should be noted that 
the incorporation of all the 4p4h effects in the present formalism requires a full four-phonon calculation. For this 
reason, a complete resolution in terms of a microscopic description of the spectrum may only be partially successful. 

We note that there has been a tremendous progress in recent years in the microscopic description of p shell nuclei 
using Green's Function Monte Carlo and no-core shell model methods |2^ [s^, |^. A possible application of the 
no-core shell model to ^^O would properly include such 4p4h effects. However the description of spectroscopic factors 
would still require the construction of effective operators to include the effects of short-range correlations on these 
quantities whereas these are automatically included in the SCGF method. 

The paper is organized as follows: Sec. ^1 introduces the formalism to account for two-phonon coupling in the 
calculation of the ph propagator. The approach employed here is based on a formalism first introduced by Baym 
and Kadanoff for the description of response functions in a many-body system at finite temperature [3^ l33l l34| . 
This framework provides a procedure to construct the effective interaction of the ph Bethe-Salpeter equation that 
generalizes the (D)RPA approach. This method is described in Sec. Ill Al and the resulting equations in Sec. IIIBI 
More technical details are left to the Appendix. Sec. IIIII describes the results for the spectrum within the current 
approximation scheme. Sec. llVl is devoted to a study of convergence properties related to the number of two-phonon 
configurations included, and the role of time-inversion diagrams. In Sec. llVXl we discuss the possible appearance of 
instabilities, in particular for the 0+ state, that was also observed in Ref. j9j]. Conclusions are drawn in Sec. IVI 
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II. EXTENSION OF ph(D)RPA FORMALISM. 

In this work the central quantity of interest is the two-time polarization propagator, whose Lehmann [ssf represen- 
tation is given by 



In Eq. 1^ , the poles and residues contain the information on the response and excitation energies of the system with 
A particles in terms of the quantities 

= - , (2) 

where and l^*^) are the exact energies and eigenstates of the A-particle system, the subscript refers to the 
ground state and cj, (cq) is the creation (annihilation) operator of a particle in the state a. For a clearer discussion 
of the formalism employed in this work, it is useful to first give a description of the relevant contributions in terms of 
Feynman diagrams. 



A. Diagrammatic contributions 

The exact ph propagator |^ is a solution of the Bethe-Salpeter equation, depicted in Fig.|2l|3^|33- This equation 
can be written schematically as 

n = n^' + n-^' k^'p*"'^ n (3) 

where II-'' represents the free propagation of a quasiparticle and a quasihole in the nuclear medium and the ph kernel 
Xiph) jg^ general, a four-time quantity. According to the Baym-Kadanoff procedure, a solution for 11 is obtained 
by first generating a self-consistent solution of the sp propagator using a given choice of the self-energy. From the 
functional derivative of this (self-consistent) self-energy with respect to the corresponding sp propagator, one then 
obtains the irreducible ph interaction K^p^^ that generates the corresponding conserving approximation for the ph 
propagator (when used in the Bethe-Salpeter equation). The standard RPA approach is derived by applying this 
procedure to the the Hartree-Fock (HF) sp propagator and self-energy. This corresponds to approximating K^p^^ 
with the bare interaction V and employing bare (HF) sp propagators as external lines. Eq. ||2J) then generates the 
RPA series of ring diagrams shown in Fig. |2| 

The extension of the RPA formalism, proposed and implemented by Brand et al. [THl ITfil | , was suggested by the 
observed fragmentation of the sp strength and the necessity to go beyond a lowest-order self-energy for a commensurate 
theoretical description. There, the Baym-Kadanoff procedure is applied to a second-order approximation for the self- 
energy. This yields contributions to the kernel K'^p^^ that automatically include all the terms which couple the ph 
states to the 2p2h ones, in accordance with the Pauli principle for the latter states. Thus, this formalism takes into 
account the mixing with 2p2h configurations in the construction of the ph propagator (the diagram of Fig. 2Ji gives 
an example). In the work of Brand et al. only a singic-polc approximation to the self-consistent propagators was 
employed. It is the aim of the present work to take the effects of the sp fragmentation more completely into account and 
therefore a fully dressed sp propagator must be used. The one employed in the present work was obtained in Ref. 
by means of a Faddeev expansion for the nuclear self-energy. When one applies the Baym-Kadanoff prescription to 
the latter self-energy, a large set of contributions to K^p^^ is generated |33|. According to Ref jl^l, one may expect 
that the most important of these terms involve the couplings to two ph phonons as depicted in Fig. ^ It is not 
difficult to see that different diagrams, similar to those in Fig. can be obtained through Pauli exchange of the 
phonon's external lines. In total there are sixteen such possible contributions, corresponding to all the possibilities of 
connecting two ph phonons to a ph state by means of a single interaction, both in the upper and lower part of the 
diagram. It is this approximation to the irreducible interaction that will be pursued in the present work. 

An additional ingredient entering the 2p2h Extended RPA (ERPA) of Refs. 0,0, "3^ requires further discussion. 
This involves diagrams similar to the one in Fig. ^ that are obtained from the 2p2h contribution by inverting the 
sense of propagation of either the incoming or the outgoing ph pair. These diagrams involve higher excitations (at least 
3p3h, when combined in the expansion of Fig. |2J) and are expected to give rather small contributions. Nevertheless, 
they represent corrections to the terms of the ph interactions that control the RPA correlations and they also add 
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Pauli corrections to the RPA expansion of Fig. at the 3p3h level. In Ref. Tej it was found that they play a role in 
stabilizing some particular solutions. 

The present calculation includes both the direct two-phonon contributions of the type depicted in Fig. and 
the diagrams similar to the one in Fig. Obviously, all the 2p2h ERPA contributions are incorporated in this 
approach. Moreover, full two-phonon configurations are accouncd for and a SCGF ap pro ach is a ppli ed. Thus, the 
present formalism is an extension of (and goes well beyond) the calculations of Refs. [THl Isfll l4f]l| . However, to 
avoid complications in the notation, we will still refer to it as "Extended RPA" , or as "two-phonon ERPA" whenever 
confusion may arise. 

The Baym-Kadanoff procedure also generates other two-phonon contributions, for example those coupling a pp and 
a hh phonon. These could be mixed with the presently considered configurations by means of an all order expansion 
of the Faddeev-Yakubovsky type. Given our present knowledge, such a massive resummation of diagrams does not 
appear to be relevant for the understanding of the spectrum of ^^O. Such a study is in any case beyond the scope of 
the present paper. 



B. Two-phonon contributions to the ph propagator 

The usual DRPA equations are obtained from Eq. Q by choosing Kj^^ l^g = VaS.js-y, as mentioned above, and by 
keeping dressed propagators as external lines. In this way, one is left with only two-time quantities and, after Fourier 
transformation, the DRPA equation becomes 

n(u;)a/3,7<5 = ^Lp.nsi^) + n^/3,Mp(^) ^^'^.P'' n^a.nsi^) , (4) 

where all the indices and summations arc shown explicitly. In Eq. the free polarization propagator (ui) is also 
a two-time quantity, with the following Lchmann representation 

where = i^n^^l'^al'^o) (3^a = (^fe "^kal^cf)) spectroscopic amplitudes for the excited states of a system 

withA-|-l {A— I) particles and the associated poles e+ = E^^^ — {e^ = Eq — E^^^) correspond to the excitation 
energies with respect to the A-body ground state. The indices n and k label the eigenstates of the systems with A±l 
particles and enumerate the fragments associated with the one-particle and one-hole excitations, respectively. When 
a dressed propagator is used as input, its one-body overlap functions and quasiparticle energies already contain 
information about the coupling of sp motion to 2plh, 2hlp and more complex configurations. As a consequence, 
contributions beyond the Iplh case are already included in n'''(a;). 

Methods based on an RPA- like expansion produce an infinite series of diagrams in which the direction of propagation 
can be reversed from backward to forward and vice versa. In the case of standard (D)RPA, the interaction kernel is 
simply given the potential V and it is the same for every contribution to the diagrammatic expansion. Therefore the 
usual (D)RPA equation can be written in compact form, as in Eq. |0J. This is no longer true when one aims to include 
additional contributions and at the same time insists on working with two-time quantities. The two-phonon diagrams 
of Figs, n and ^ have different analytical expressions. As a consequence one first needs to separate all the four 
possible time directions of the kernel if (p'') -forward to forward, backward to backward, and the two time-inversion 
cases- before including the relevant diagrams in the Bethe-Salpeter equation 0. This can be achieved by splitting 
the free ph propagator |(5j into its forward- and backward-going parts, denoted by ^ and ^ respectively, 

Wiu) — > >iu;) + <H . (6) 

By performing this substitution in Eq. (j2Jl one obtains a similar separation for the complete ph propagator 

u{uj) ^n>(c^) + n<(w) , (7) 

where ^ and ^ now refer to the sense of propagation of the final lines only. Suppressing the indices and summations 
one obtains 

Il>{u;) = n-^ >{u;) + >{iu) K^P''^ n(w) , 

Il<{u;) = <(w) + <{lu) K^p''^ II{lu) . (8) 
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The last step consists in substituting Eq. ^ into (jSj) and approximating each component of the ph kernel K'^p^^ 
with the sum of the bare interaction and the corresponding two-phonon contributions. The result corresponds to the 
ERPA equations given by 

n>(w) = n-^>(w) + nf>{uj) {{v + w>{u))\i>{uj) + {v + i7>'<)n<(w)} 

n<(cj) = n/<(cj) + <(uj) {{v + H<'>)n>(uj) + {v + w<{uj))u<{uj)] . (9) 

In Eq. W^{Ld) represents the contribution of all the sixteen two-phonon diagrams, Fig.^ in the forward direction. 
VF^(w) corresponds to the contributions connecting the backward-going terms. Analogously H^'^ represent the sum 
of the diagrams in which a backward-going hp state is inverted in a forward-going ph one (Fig. ^ illustrates such a 
case) while is the time-reversed contribution. 

The practical implementation of Eqs. © requires additional manipulation to treat the freely propagating lines in 
Figs, n and This situation is completely analogous to the one already discussed in Ref. 28, 41] for the 2plh 
expansion of the self-energy. Specific details for the ERPA equation Q are given in the sections in the Appendix. 
The next section reports on the application of this formalism to ^^O. 

III. RESULTS 

As in the work of Ref. the ERPA equations where solved within a model space consisting of a finite set of 
harmonic oscillator states. All the first four major shells (from Is to 2plf) plus the 159/2 where included to account 
for the sp orbitals that are most relevant for low-lying excitations. The harmonic oscillator parameter was chosen to 
be 6=1.76 fm. As a consequence of the truncation of the model space a Brueckner G-matrix was used as a microscopic 
effective interaction. This G-matrix was derived from a Bonn-C potential and computed according to Ref. ^_43j . 

The contributions of two-phonon states were first studied by solving the ERPA equation @ with an independent- 
particle model (IPM) propagator. This propagator was constructed from a Slater determinant composed of the lowest 
occupied harmonic oscillator (h.o.) orbitals in the model space. Where the effects of nuclear fragmentation were 
included, the calculations employed the fully dressed sp propagator computed in Ref. 0. This propagator contains 
no more than two main quasihole (quasiparticle) fragments for each state in the p (sd) shell. For these orbits, the sp 
propagator can therefore be well represented by means of one or two principal fragments and an effective pole that 
accounts for strength far from the Fermi energy. For quasiparticle states associated with the pf shell, instead, the 
fragmentation is more substantial M and the self-consistent propagator was approximated by including the most 
important poles (up to 4 main fragments for the P3/2) plus two effective ones that gather all the background strength. 
We have checked that the results are not sensitive to the details of the fragmentation of the If orbitals. For the 2p 
orbitals, we found that the strength distribution is sufficiently represented by the adopted prescription and further 
discuss the sensitivity to the results in Sec. IIV Al 

In both cases (IMP and dressed inputs), the standard (D)RPA equation was solved first, in order to generate the 
ph phonons that entered the final two-phonon ERPA calculation. When coupling different phonons, only the lowest 
few states of ^^O have the right quantum numbers, and energies low enough, to generate two-phonon contributions 
with unperturbed energy below 20 MeV. These are the onry states relevant for our purposes 0. In practice, 
the calculations were performed by including, in each channel, all two-phonon configurations up to 30 MeV. The 
stability tests reported in Sec. IIV 01 demonstrate that, for the case of a dressed input propagator, higher two-phonon 
contributions do not change the results for the low-lying states appreciably. 

A. Results for the particle-hole propagator 

The results obtained from the h.o. input propagator are displayed in Fig. [S] The isoscalar eigenvalues obtained 
for energies below 15 MeV are displayed, for both the standard RPA and for the ERPA calculation. The ERPA 
calculation produces lower energies for some of the states that were already obtained in RPA. Both 3^ and 0+ 
shifted down to about 0.5 MeV above the experimental energy. As discussed further in Sec. IIV CI this is to be 
considered rather fortuitous and one must remember that the dressing of the sp propagator, which has the effect to 
raise these eigenvalues, has not been taken into account yet. Still, this result indicates that correlations between ph 
and two-phonon states can be substantial and go in the right direction in explaining the experimental results. 

It is also worth considering the total ph spectral strength of these states obtained by summing the corresponding 
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amplitudes, Eq. |5J, as follows 



= E p5 • (10) 

a/3 



Results for both the excitation energy and of the principal levels in Fig. are given in Table HI For the ERPA 
case, also the fraction of quasiparticle-quasihole and two-phonon configuration that appear in the wave function are 
shown. Note that is substantially bigger than one for the 3~ and 0+ states and that these values increase further 
for the ERPA results. This signals an increase of the collective character of these solutions which may lead to an 
instability of the RPA equations for interactions that are even more attractive. The two-phonon ERPA approach 
generates a triplet of states at about 14 MeV, with quantum numbers 0+, 2+ and 4+. However, these levels are 
almost exclusively composed of two-phonon contributions and contain only small admixtures of ph states, resulting 
in a small ph spectral strength . The quantum numbers and energies of these states indeed correspond closely 
to those obtained by coupling two 3^ RPA phonons, each at 7.14 MeV. We note that a similar triplet is found 
experimentally at 12.05, 11.52 and 11.10 MeV, which also corresponds to twice the experimental energy of the first 
3~ phonon. The first experimental 2"*" is found at lower energy and its spectral strength is known to have relevant 
ph components Thus it cannot be identified with any of the above two-phonon contribution. For all the lowest 
states that are not already reproduced by standard RPA a very small ph spectral strength has been found, due to a 
general lack of mixing between the ph and the two-phonon configurations. Of interest is also the 2"*" that represents 
the giant quadrupole resonance at 20.7 MeV. In this case RPA and ERPA give 22.9 and 23.3 MeV for the main peak 
but with a lower Zn„ in the second case. For this state, part of the ph strength (about 20%) is shifted to two-phonon 
configurations representing the expected fragmentation of the giant resonance. 

Fig. and Table HT) show the analogous results when dressed sp propagators are used as input. In general, the 
main effect of fragmentation is to screen the nuclear interaction as a consequence of the quenching of spectroscopic 
factors for the input sp propagator. For the DRPA case, this results in increasing the lowest 3~ and 1~ solutions 
by ~2 MeV. More substantial is the effect on the lowest 0+ state with a predominantly ph character which rises to 
about 17 MeV, confirming the sensitivity of this state to details of the fragmentation and the strength of the nuclear 
interaction. Unlike the IPM case, we have chosen to solve the two-phonon ERPA by first shifting the lowest solution 
for both the 0^, 3~ and 1~ states down to their relative experimental energies. This has the advantage of lowering 
the most important two-phonon configurations and allows to investigate their interplay with the ph ones. 

As can be seen in Fig.|Hl the ERPA equations still generate a triplet of 0"'', 2+ and 4+ states at twice the energy of 
the first 3~ phonon. Due to the screening of the ph interaction, these states mix very little with the ph configurations 
yielding an almost degenerate triplet. Table Hll also gives a comparison between the total two-phonon content of the 
states and the individual contributions of the most important configurations. This decomposition demonstrates that 
this triplet is formed by pure 3~ ® 3~ states. This observation is in accordance with the 2p2h character of these 
states 23, 45]. Since the overall energy of these states approximately agrees with experiment, one can expect that this 
calculation correctly represents the bulk properties of their wavef unctions. Nevertheless, it is clear that an additional 
interaction is needed in order to split this triplet as observed experimentally. It should be noted that no solution that 
can be identified with the 2+ and 4+ levels at 9.8 and 10.4 MeV. The 2^^ level has been interpreted in terms of an 
a particle rotating around an excited ^^C core "23^ and therefore involves correlations that may go well beyond the 
present calculation. With regard to the 2+ strength around 20 MeV, we note that the sp fragmentation included in 
the DRPA equation already generates a distribution of 2+ strength. Table ^ contains a few of these solutions and 
shows that also for the dressed case a sizable mixing of ph and two-phonon contributions is obtained at this energy. 
This mixing generates the spreading of ph strength over different solutions, consistent with the finite width of such 
a resonant state. It should be noted, however, that a discrete basis is used here and therefore it is not possible to 
properly describe the continuous strength distribution of a giant resonance, for which a continuum or a complex basis 
should be used [46l |. 

Also interesting are the results for the other isoscalar 0''' states. The lowest solution with a predominant ph 
character, that was found in DRPA at ~17 MeV, is obtained at a similar energy but is now characterized by a partial 
mixing with two-phonon configurations. The ph and two-phonon contents of this state given in Table HTl show that this 
is the result of mixing with the the lowest solution, which ends up at ~11 MeV. The latter is predominately a two- 
phonon state. It is also seen that in both cases the relevant configuration comes from the coupling of two 0^ phonons. 
Configurations involving two 3~'s make smaller contributions while the other modes not reported in the table were 
negligible. We note that the wavefunctions for these states contain several relevant ph configurations, obtained from 
different quasiparticle fragments in the pf {sd) shells combined with quasihole fragments of the p (s) shells. Therefore, 
the situation is more complicated than the simple picture of only two levels interacting with each other. We observe 
that the shell-model calculations for the first excited state of ^^O give very small contributions from both Qhcu and 
2huj configurations and a strong population of 4huj states . On the other hand, inelastic electron scattering 
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experiments [23j clearly excite this state identifying its partial ph character even though the state is dominated by 
4p4h components. From the point of view of the SCGF approach, the one-body response is completely described by 
the (dressed) ph propagator, Eq. (Q, and therefore the total one-body strength must be represented by ifTIUl . This 
observation points to the necessity of a stronger mixing between ph and two-phonon configurations than obtained in 
the present calculations. Still, the present results indicate that the two-phonon configuration obtained by coupling 
the two lowest 0+ phonons plays a role in the structure of the first 0+ state itself. This means that, starting from 
the dressed ph admixture contained in it, a fully self-consistent calculation would also generate some contributions of 
2p2h, 4p4h, and beyond. The calculations of Ref. 25] have shown that the bulk of the 4p4h contributions to the first 
excited state of ^^O may come from the coupling of four different phonons with negative parity (3~ and However, 
the self-consistent role of coupling positive parity states was not conidered in that work. Both this effect, the RPA 
correlations and the inclusion of the nuclear fragmentation allow for the -at least partial- inclusion of configurations 
beyond the 2p2h case even when no more than two-phonon coupling is considered, as in this paper. A study of the 
importance of three- and four-phonon configurations within this approach is beyond the scope of the present work. 

The principal negative parity isoscalar states are still shifted down in energy by ERPA over DRPA, as expected, but 
the improvement is less than 0.5 MeV. The low-lying 3~ and 1~ levels remain substantially above the experimental 
energy at 9.23 and 10.90 MeV, respectively, and more correlations will be needed in order to lower their energy. We 
note that these wavefunctions contain a two-phonon admixture obtained by coupling the low-lying 0^ excitation to 
either the 3~ or 1~ phonons. According to the above interpretation of the first 0"*" excited state this corresponds 
to the inclusion of 3p3h and beyond. At the same time two additional negative parity solutions, with two-phonon 
character, are found at higher energy in accordance with experiment. 

The ERPA results for an h.o. input give appropriate corrections to several energy levels. However, such corrections 
tend to become negligible in the successive calculations employing dressed propagators. Presumably this happens 
because a sizable part of the correlations, which in the latter case are introduced by ERPA, are already accounted 
for by the dressing of the sp propagators. In general, the inclusion of fragmentation has also the effect of screening 
the interaction between ph and two-phonon configurations. However, the results of Table HTl show that a considerable 
mixing between the two can be obtained if all the relevant configurations are sufficiently low in energy, although 
no substantial mixing between different two-phonons states is seen. According to these considerations, the biggest 
deficiency of the present approach is probably the lack of an interaction within the multi-phonon space. Refs. |25ll26l| 
have shown that such correlations are due to pp and hh interactions, neglected here, and that their inclusion has 
important effects on the final result for the spectrum. These terms may considerably lower most of the eigenstates 
of the ERPA equations. However, a second type of correlation between multi-phonon states, that was not accounted 
for in Refs. [25ll26| . is the Pauli exchange between the fermion lines included in different phonon propagators. These 
will probably generate the opposite effect, screening part of the pp and hh correlations. 

Although the results obtained with a dressed input propagator for the lowest excited states are not completely 
successful, it must be recognized that the present ERPA calculation does represent a step forward. Indeed a proper 
description of two-phonon contributions, at the present level, had not been accomplished before. These states not 
only account for configurations at the 2p2h and higher level but appear to have some relevance for the description 
of states at higher energy. In addition, excitations are generated with quantum numbers that can not be obtained 
in the standard (D)RPA approach at these low energies. Given the results of the present work, it appears that a 
calculation including both pp and hh correlations, a more complete treatment of Pauli effects and the most relevant 
configurations up to the four-phonon states can be achieved with the presnt-day computers. Such a calculation is 
planned for the future. 

IV. STABILITY OF THE RESULTS AND TRUNCATION OF THE MODEL SPACE 

A. RPA instability 

The RPA approach may break down for strongly attractive interactions, by genereating an excess of collectivity and 
lowering the first excited state below the ground-state energy. The inclusion of fragmentation of the sp strength can 
affect this behavior in several ways. First, both the screening of the nuclear interaction and the splitting of the particle 
and hole spectral strengths over different fragments tends to stabilize the DRPA equations. Second, correlations lower 
most of the sp energies in comparison with the IPM. This results in smaller unperturbed ph energies and therefore 
has the opposite effect of pushing the DRPA approach toward instability. Third, the instability can be sensitive to 
the details of the spectroscopic amplitudes. This last point can be best illustrated by looking at the sp strength with 
j'" — 3/2^, to which the present calculation is most sensitive. For the model space employed in this work a given 
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particle overlap function is a superposition of two h.o. eigenstates, 

^;3/.W = -^I'k/. ± ^^2^3/. ^ (11) 

where we have stressed the fact the two components can sum either constructively or destructively. The fragmentation 
introduces a non-zero X^^^^^ component even for particle states that in the IPM would be described only by the 2p^/2 
subshell. By considering Eqs. (|A14|) one can see that the ph interaction for the DRPA approach can be substantially 
changed by this new component by virtue of the strong matrix elements of the effective interaction between 
states. In particular, this also applies to the off-diagonal terms that drive the RPA-like correlations. It is also worth 
nothing that the sign of this correction changes accordingly to the relative sign of the components in Eq. Hll|l . therefore 
pushing the states either towards stability or instability, respectively. 

We have observed that choosing a wrong sign for one or more of the particle fragments has drastic effects 
for the solutions of the isoscalar 0+ channel, with the choice of constructive interference in Eq. leading toward 
instability of the DRPA approach. For the input propagator employed in this work, all of the background fragments 
have a similar wavefunction and a negative or positive interference sign according to whether their sp energy is lower 
or higher than 20 MeV, respectively [s^. Therefore it is natural to collect such background distributions in two 
different poles, one for each region of missing energy. A similar behavior applies also to the other quasiparticle states 
in the 2p and sd shells. We have checked that the results of Sec. IIIII are stable with respect to the number of main 
fragments and effective poles included in each subshell, therefore no instability affects the present work. However, 
in applications in which the number of poles has to be strongly reduced it may not be possible to chose effective 
poles of the sp propagator that represent the details of the spectral distribution sufficiently accurately. In such a 
case an instability in the DRPA approach may arise. This situation occurred with a previous DRPA calculation 
of the isoscalar 0+ channel, reported in Ref. (9j, in which a smaller number of fragments was included. Thanks to 
improvements in the computer code that solves the DRPA problem, it is now possible to consider a sufficiently large 
number of poles in the sp propagator and this problem has been overcome. 

B. Time-inversion diagrams 

The contribution of two-phonon terms to the time inversion diagrams of Fig. 0] can, in principle, generate other 
correction to the interaction between ph and hp states that drives the RPA correlations. Nevertheless, the effect is 
very small due to the large energy denominators that appear in Eqs. ljA18|l and ljA19|) . We have tested their influence 
by neglecting the corresponding terms and H^'^ in Eq. 10 and report the results in Table UllI No appreciable 

change is generated illustrated by differences with respect to the results of Sec. lIIII of at most 0.1%. It is worth noting 
that the contribution of H^'"^ and H"^'^ can in principle carry information on the Pauli breaking at the level of 3p3h 
and beyond. Therefore, they may become more important for the case of three-phonon calculations. Nevertheless, 
Table Hill suggest that they are not likely to play an important role in the description of the spectrum of ^^O. 

C. Stability vs. number of two-phonon configurations 

Fig. [3 shows the results for selected solutions of the ERPA equations, obtained by employing different sets of two- 
phonon states. For any given point, only those configuration with energy ej^^ -I- e^^ < Ecut have been included in 
the calculation. For Eciit=70 MeV, about 700 to 1000 two-phonon contributions have been included, depending on 
the channel. This is one order of magnitude larger than the number of ph states that enter the calculation. Since 
no two-phonon configuration has energy lower than 10 MeV, the leftmost points in Fig. [7| correspond to the simple 
(D)RPA calculation. As the lowest few two-phonon contributions are included, the solutions for these levels show a 
sudden jump for all the lowest excited states. After this, the results obtained by using an IPM input still continue to 
exhibit a dependence on the value of the cut-off, roughly lowering states by 1 MeV every time that Ecut increases by 
20 MeV. In particular, the 0+ and 3~ solutions decrease progressively, eventually heading to RPA instability. This 
confirms that the values reported in Tabled for these states do not correspond to converged results. The situation is 
much better for the case of a dressed input propagator, for which the low-lying solutions are approximately stable for 
Ecut > 20 MeV. This confirms that corresponding higher-energy two-phonon excitations do not have strong influence 
as their effects have already been included by the dressing of the sp motion. 

In Fig. [3 we also show the results for the most relevant of the 2+ solutions that represent the resonance at about 
20 MeV. Here the variation is more appreciable, in particular for values of the cut-off Ecut comparable with the 
energy of the state itself. Also here a more stable behavior is obtained for the case of a dressed input propagator. We 
conclude that when the nuclear fragmentation is accounted for the choice of Ecut = 30 MeV adopted in Sec. IIIII is a 
adequate for considering the low-lying spectrum. 
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D. Size of the model space 

As a second possible source of uncertainty one may consider the size of the employed model space. The choice 
of the model space employed in Sec. IIIII allows to take into account excitations up to the 2plf shells. This means 
that ph configurations up to 2huj ~ 27 MeV are included, consistently with the choice for truncating two-phonon 
states discussed above. The calculations of Ref. 0, [s^ also show that this is adequate and suggest that better 
convergence is to be expected for a calculation with a dressed input. To check this for the present calculation, we 
have augmented the model space by adding two more h.o. shells, up to p = 2n + €=5, and recomputed the G-matrix 
elements accordingly. The input sp propagator used in the calculations was obtained by including an extra particle 
pole for every subshell added to both the IPM's Slater determinant and the fragmented dressed propagator. The sp 
energies for the new shells were chosen by solving the Brueckner-Hartree-Fock equation within the new model space, 
in an approach analogous to the one of Ref. @ . 

The results are compared to the ones obtained for the smaller space in Table Hill For the IPM, a sizable variation 
of the 3~ excitation energy is found. Also, the high values of the ph strength for both the 0"*" and 3~ solutions 
indicate that these states are approaching instability. As a consequence of the different 3~ excitation energy, which 
also enters the two-phonon calculation, most of the ERPA solutions are shifted down in energy. The main conclusions 
of Sec. IIIII however, remain unchanged. When dressed sp propagators are employed, no dramatic change occurs for 
the low-lying solutions and only some states with dominant ph character are shifted by no more that 1 MeV. This 
suggests that the main conclusions of Sec. IIIII are also not affected by the truncation of the model space, whence the 
bulk of the missing correlations has to be looked for in an extension of the diagrammatic expansion. We note, however, 
that a more attractive G-matrix is expected to be obtained if the continuum outside the model space were modified 
by accounting for self-energy corrections. This feature may lead to a further improvement of the present results. In 
addition, one may explore the inclusion of effective three-body force as due to propagation outside the model space. 
Such effective forces are critical in obtaining the correct binding energy in no-core shell model calculations reported 
in Ref. [sj. 

V. CONCLUSIONS 

The dressed RPA equations have been extended to account for the coupling of two ph phonons in forming the excited 
states of a many-body system. The coupling among ph and multi-phonon configurations is conceptually similar to the 
Interacting Boson Model of Ref. ^25] but limited, in this application, to only two-phonon admixtures. Nevertheless, 
the present approach has the added advantage of taking into account both the effects of nuclear fragmentation and 
the RPA-like correlations, as generated by two-phonon fluctuations in the ground state. 

The resulting Extended RPA formalism has been applied to study the excitation spectrum of ^^O. The results 
suggests a sizable mixing of ph and two-phonon configurations for the case of the 2+ quadrupole resonance at 20.7 MeV. 
Other solutions, carrying quantum numbers that cannot be generated by the simple DRPA, are obtained with this 
method, among which a few isoscalar negative parity states and a triplet with 7^=0"*", 2+ and 4+ near 12 MeV, in 
accord with experiment. In particular, the states in the triplet were seen to be almost pure two-phonon configurations 
obtained by coupling two 3~ phonons. This feature confirms the 2p2h character of these states. The results are less 
satisfying for the low-lying positive parity states, which are known to require a proper description of 4p4h excitations. 
The present approach predicts the lowest solution for the isoscalar 0+ channel to be at 11.3 MeV, considerably above 
the first experimental excited state. A sizable component of this solution is seen to be generated by the coupling of two 
of the lowest 0^ phonon themselves. This suggests that important contribution of 4p4h excitation may be included 
already at the two-phonon level in a self-consistent fashion. Other correlations beyond the 2p2h level are partially 
included in the present calculation through the dressing of the propagators and the RPA approach. As expected, 
though, the present implementation is not sufficient to obtain a complete description of the low-lying excitation 
spectrum of ^^0. The results obtained in this work show that an interaction between multi-phonon configuration 
is still missing. This can be achieved by including pp and hh correlation between different ph phonons. Also the 
inclusion of up to four-phonon states is expected to be relevant for this system. 

Finally, the stability of the present results has been tested with respect to the truncation of the model space and to 
the number of two-phonon configurations accounted for. It was found that the effects of nuclear fragmentation acts to 
'renormalize' the sp propagator, making the solutions of the ERPA fairly independent of higher energy configurations. 
As discussed in Sec. IIVI this feature generates stable solutions with respect to the size of the model space and the 
corresponding G-matrix interaction used in this work. 

The four-phonon calculation of Refs. |2^ |2^ give a very good description of the excitation spectrum of ^^O. 
Those findings and the effects of fragmentation discussed in this paper suggest that the present ERPA formalism 
can be suitably extended to generate a satisfactory description of this nucleus, within the framework of SCGF. Such 
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extensions are presently under consideration. Such a calculation may also significantly reduce the discrepancy between 
the measured and the theoretical sp spectral function for this nucleus 0. 
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APPENDIX A: EXTENDED DRPA EQUATION WITH TWO-PHONON CONTRIBUTIONS 

When reducing the kernel of the Bethe-Salpeter Eq. © to a two-time quantity, one has to deal with the fact that 
some lines in the diagrammatic expansion continue to propagate unperturbed while some interaction occurs between 
other particles. This is also the case for the diagrams of Figs.nandQl This situation can be overcome by redefining the 
objects that appear in Eq. © and promoting the quantum numbers {n, k}, -that label particle and hole fragments- 
to external indices. The usual form of the ph propagator, Eq. is obtained only as a last step of the calculation by 
performing the sum over the {n, k} indices |28ll39| . 

The separation of the propagators n'^(a;) and n(a;) into forward and backward components, outlined in Sec. IIIBI 
is a natural consequence of adopting this prescription. Indeed particle and hole external lines turn into each other 
by time inversion and become quantities that depend on different quantum numbers, i.e. the fragmentation indices 
{n, k}. For this reason, Eqs. ((BJ and ((TJ are only formal relations and the arrow can be substituted by an equal sign 
only before summing over all the particle and hole fragments. The free polarization propagator @ naturally splits in 
two components that are purely forward and backward-going. The relevant Lehmann representations are 



Gt — _ G , (Al) 



nf < r ^ - _ A A ->^a ^13 [y-f Xg ) 

^ ~r \t-nfi ^ka) I 

^ G* , (A2) 

where no summation is implied and a shorter notation for the unperturbed poles and residues has been introduced. 
The quantities D and G can be thought as matrices whose elements contain all the unperturbed poles and residues, 
respectively. Note that in this case D is diagonal and depends on the fragmentation indices {n, k} only, while G is 
rectangular because its column indices depend on the model space orbitals {a} as well. 

The separation of the complete propagator is a little more complicated. The splitting of Eqs. {7)) and |(SJ) involves 
the time direction of the outgoing lines only. Since the RPA series contains contributions that can invert several times 
the sense of propagation of the ph diagrams, both forward and backward poles can appear in each component 

{ 7> n A -p-n 57< n ( yn \ 

iia„./3fc,,76M - u^-el + i^ ^ uj + el-ir^ ' ^^^^ 

/ y< n \ yn y> n ( yn \ 

n^O " ' n^a " ' 

In Eqs. I|A3|I and (|A4(1 the spectroscopic amplitude splits in two contributions and These appears unchanged 
in both equations due to the time-inversion symmetries obeyed by the ph states. In terms of these definitions, and 
applying the summation prescription, Eq. iQ can be exactly formulated as follow 



yn _ 



rii ,k2 
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"1 ,^2 



(A6) 



Finally, the contributions of two-phonon diagrams in the forward and backward direction can be expressed in 
Lehmann representation as well 



(-) = E 



K 



>t 



1 

iu — E + irj 



(A7) 



akaf3ns,^k-y5nsy^ ) / j i /' tt i tt ^ 



= A-<t ^—K< , (A8) 

uj — bj — w] 



where eJJ and are the energies of the intermediate ph phonons. 



ERPA Matrix 



The eigenvalue equation for the ERPA is obtained in the usual way, by substituting Eqs. (|A1|I to l|A4|l into Eq. 
and then extracting the residues of the ph poles e^. The result is an eigenvalue equation in terms of the vectors 
and Z"^ . 

To linearize the problem, we introduce the following components 
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where X^^'^ and y*-^-* represent the ph amplitudes that appear in the standard (D)RPA equations [la. |47| and AT^^^ 
and 

are the analogous two-phonon amplitudes introduced by the ERPA approach. The components X^^^ and 

yd) 

are related to Z-^ and respectively by 
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Eqs. © can be put in the form of a linear eigenvalue equation 
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where the matrix M is defined as 
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The off-diagonal 2x2 blocks in Eq. ljA12|l describe diagrams in which the time direction of propagation is inverted. 
In the present case the only non vanishing elements are the ones that involve the inversion of a single ph state into a 
hp one, or vice versa. These correspond to the sum of the first-order term V, which represents the kernel of the bare 
RPA, and the more complex diagrams of Fig. ^ . Blank spaces would in principle allow to include more complicated 
contributions, that involve time inversion of two-phonon diagrams into ph or hp configurations. These contributions 
are not expected to play a relevant role for the present problem. It must be noted that if the terms involving the 
matrix G are discarded in Eq. ljA12|) . the components X^^^ and F^^^ decouple from X'^^ and F^^). In this case, 
Eqs. © reduce to the ph-DRPA one Q and the matrix (jA12p would take the form of the standard RPA matrix 41^. 
The normalization condition is derived in the usual way, by extracting the contribution of order zero of the expansion 
around a given pole and by employing the conjugate equation. One eventually obtains |16| | 

J2 (^^''^ - Y'^'^^ yW) = 1 , (A13) 

4=1,2 

where the inner product of the vectors X*^*^ and F^^'^ is implied. 



2. Matrix Elements for ph ERPA 



In the following we give the explicit expression for the matrix elements of Eq. ljA12p . The contributions originating 
from the standard (D)RPA equation are (here and below, summations over repeated greek indices are understood) 



(GVG^) , , 
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112 ,n3ki 
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with the corresponding unperturbed ph energies 

The interaction between two-phonon intermediate states and the ph ones is given by 

(^^'>Gt)„.„,„^,.^ = \ {[(x(i)„:ix(\l)* x-^y-^^v,.,.y^,^{yl^-) 



(A15) 



(A16) 



where ri^ and are the quantum numbers of the two phonons that form the intermediate state. The quantities 
X(^)„"j,2 are the forward-going amplitudes of the intermediate phonons (|A9(I or l|A10|l and are obtained from the 
previous solution of the polarization propagator. Note that K^G'^ is symmetric under the exchange of the two indices 
n1 and n^, as required by the boson character of the ph phonons. The factor ^ assures that no double counting 
happens when only two free propagators 11^ (w) are coupled. Due to the time- inversion properties of npnh states, the 
corresponding contribution for backward-going propagation is simply related to above one: 



{K< (G*)^) 



{(K>G^) „ , , }* 
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The analytical expression for the contribution of the time- inversion diagrams is a little more complicated. First, 
we introduce the following quantity that corresponds to the generation of a two-phonon state by a fluctuation of the 
nuclear mean field 
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This is also symmetric under the exchange of the indices and n^. Then, the matrix elements containing the 
time-inversion diagrams, H-''^ and H^'-^, can be written as follows (a summation over and is also implied) 

+ [- («! -s-^ "4) - {k2 ^3) + {ni n4, fc2 kz)] 



+ [ni 714,^2 ^3] 



Finally, the two-phonon unperturbed energies are given by 
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FIG. 1: Examples of contributions involving the coupling of two independent ph phonons. In total, there are sixteen possible 
diagrams of this type, obtained by considering all the possible couplings to a ph state. The two-phonon ERPA equations Q 
sum all of these contributions in terms of dressed sp propagators. 
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FIG. 3: Diagrammatic expansion of the standard RPA equation. An explicit time direction is assumed for the diagrams of this 
figure. 
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FIG. 4: Example of direct, a), and time inversion, b), diagrams that appear in the standard ERPA expansion. Both diagrams 
a) and b) come from the same four-time screening diagram. The last picture in part c) shows the corresponding two-phonon 
extension of the time-inversion contribution. Note that the diagram b) generates Pauli exchange corrections to the last diagram 
shown in Fig.|3 
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TABLE I: Excitation energy and total spectral strengths obtained for the principal solutions of the RPA and ERPA equations. 
For the ERPA case the total fraction of ph and two-phonon contributions are also shown. An IPM input sp propagator was 
used to generate these results. 
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spectral strength Zn„ lower than 10%. The first column reports the experimental results taken from Ref. j44|. 
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FIG. 6: Results for the DRPA and the two-phonon ERPA propagator of ^^O with a dressed input propagator from Ref "9^, 
middle and last column respectively. In solving the ERPA equation, the lowest 3~, 1~ and 0"'' levels of the DRPA propagator 
where shifted to their experimental energies. All other DRPA solutions were left unchanged. The excited states indicated by 
dashed lines are those for which the (E)RPA equation predicts a total spectral strength Zn^ lower than 10%. The first column 
reports the experimental results f44j . 
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FIG. 7: Dependence of the ERPA solutions on the number of two-phonon states considered. For any given point, all the 
configuration with energy ej^^ ^ Ecut have been included in the calculation. Solid (dashed) lines refer to results obtained 

from a dressed (IPM) input propagator. 



r = dressed/DRPA dressed/ERPA (0+)^ (3")^ (0+,3") (0+,l") 

r el el Z^^ ph{%) 2n(%) (%) (%) (%) (%) 

2+ 23.77 0.468 23.52 0.123 26 74 

2+ 22.96 0.341 78 22 

2+ 20.59 0.269 20.42 0.255 98 2 

1" 13.37 0.148 21 79 79 

3" 12.35 0.113 16 84 84 

0+ 12.15 0.001 1 99 3 96 

4+ 12.14 0.007 1 99 99 

2+ 12.12 0.008 1 99 98 

0+ 16.62 0.717 17.21 0.633 88 12 10 0.5 

0+ 11.28 0.092 12 88 85 2 

1" 11.19 0.720 10.90 0.680 94.1 5.9 5.8 

3" 9.50 0.762 9.23 0.735 95.9 4.1 4.0 



TABLE II: Excitation energy and total spectral strengths obtained for the principal solutions of DRPA and two-phonon 
ERPA equations. A dressed sp propagator was employed. The total contribution of ph and two-phonon states of the ERPA 
solutions are shown. For states below 15 MeV, the columns on the right side give the individual contributions of all the relevant 
two-phonon contributions. The sum of these terms for the states that are listed does not exceed 1%. 
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0.008 


12.12 


0.008 


12.05 


0.013 


4+ 


13.87 0.028 


13.87 


0.028 


11.81 
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0.092 
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9.20 1.075 
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1.024 
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10.90 
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10.89 
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10.64 


0.723 


3" 


6.55 1.269 


6.54 


1.266 


4.68 


1.581 


9.23 


0.735 


9.22 


0.735 


8.91 


0.762 


0+ 


6.52 1.820 


6.36 


1.596 


6.74 


3.055 















TABLE III: The excitation energy and total spectral strengths discussed in Sec. IIIII are compared to the solutions of the 
two-phonon ERPA equations by neglecting the time-inversion diagrams of Fig. |1] Both IPM and dressed input propagators 
cases are displayed. Results obtained within a larger model space {p—2n + 1—5) are also shown. 



